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ABSTRACT 

It is shown that compact bodies project out strands of concentrated dark matter filaments henceforth simply 
called hairs. These hairs are a consequence of the fine-grained stream structure of dark matter halos, and as such 
constitute a new physical prediction of ACDM. Using both an analytical model of planetary density and nu¬ 
merical simulations utilizing the Fast Accurate Integrand Renormalization (FAIR) algorithm (a fast geodesics 
calculator described below) with realistic planetary density inputs, dark matter streams moving through a com¬ 
pact body are shown to produce hugely magnified dark matter densities along the stream velocity axis going 
through the center of the body. Typical hair density enhancements are 10 7 for Earth and 10 8 for Jupiter. The 
largest enhancements occur for particles streaming through the core of the body that mostly focus at a single 
point called the root of the hair. For the Earth, the root is located at about 10 6 km from the planetary center with 
a density enhancement of around 10 9 while for a gas giant like Jupiter, the root is located at around 10 5 km 
with a enhancement of around 10 11 . Beyond the root, the hair density precisely reflects the density layers of 
the body providing a direct probe of planetary interiors. 


1. INTRODUCTION 

As the evidence of dark matter’s existence has become over¬ 
whelming 1 , understanding its nature and interactions through 
the observation of both direct (from Earth-bound experi¬ 
ments) and indirect (excess y-ray and positron measurements) 
dark matter detection has become an increasingly vigorous 
field of inquiry. In addition, upper-limits on dark matter inter¬ 
actions are now extracted from the analysis of large samples 
of colliding clusters (Harvey et al. 2015), and constraints on 
both dark matter stream velocity dispersions and halo forma¬ 
tion can be deduced from large structures like the ’’Field of 
Streams” (Belokurov et al. 2006). In spite of all these extraor¬ 
dinary advances, little is known about the dark matter particle 
itself (whether it is a boson, Dirac or Majorana fermion), its 
local density and velocity distribution in the solar system, or 
its interactions. This missing information may remain hidden 
until a direct detection is made and their interactions can be 
studied consistently over time; a longitudinal study could be 
performed thanks to dark matter production at an accelera¬ 
tor or because a naturally occurring concentrated dark matter 
source is discovered locally that can be reliably accessed with 
instruments. Highly concentrated dark matter hairs would be 
such a source and are in fact a prediction of ACDM, a model 
that fits very well with the most recent Planck CMB analysis 
and many other astrophysical data 2 . 

For thermally produced Cold Dark Matter (CDM) of the 
type allowed by ACDM, the CDM primordial thermal ve¬ 
locity dispersion is expected to be greatly suppressed as the 
universe expands and the CDM collisionless gas cools. In 

1 In particular, the accumulation of evidence for dark matter’s existence 
spans generations (Zwicky 1933; Ostriker and Peebles 1973; Ostriker et al. 
1974; Fabricant et al. 1980) with more recent and striking observations such 
as the bullet cluster (Clowe et al. 2006). For a discussion of the challenges 
facing the Modified Newtonian Dynamics (MOND) alternative to dark mat¬ 
ter, see Dodelson (201 1). For a look at the constraints on warm dark matter, 
see Viel et al. (2013) and Anderhalden et al. (2013). For reviews of the evi¬ 
dence of dark matter, see Bertone (2010) and Massey et al. (2010); see also 
the latest Planck results for the cosmic microwave background evidence for 
cold dark matter Ade et al. (2015). 

2 See section 5 of Ref. Ade et al. (2015) and references therein. 


particular, for a WIMP with mass 100 GeV that decoupled at 
10 MeV, the velocity dispersion is about 3xl0“ 4 m/s while the 
velocity dispersion for a 10 yueV axion is 10 -8 m/s (Sikivie 
1999). As the non-linear effects of gravity become more 
prominent and the halos grow, a coarse-grained veloc¬ 
ity dispersion of the CDM will appear as they orbit the 
galaxy (Sikivie et al. 1995); the coarseness of this dispersion 
will gradually smooth out as the number of orbits increases. 
In the “Field of Streams” (Belokurov et al. 2006), the effec¬ 
tive velocity dispersion is 10 km/s providing an experimental 
upper-limit on the velocity dispersion of the fine-grained dark 
matter streams, each of which should have a tiny primordial 
velocity dispersion. The “Field of Streams” is an example 
of a tidal dark matter stream built-up from a huge number 
of fine-grained dark matter streams and distinct from them. 
A phase-space perspective sheds additional light on the pro¬ 
cesses affecting the CDM under the influence of gravity. 

At the last scattering surface, the CDM occupy a 3- 
dimensional sheet in the 6-dimensional phase space since they 
have tiny velocity dispersions. The process of galactic halo 
formation cannot tear this hypersurface filled with the colli¬ 
sionless gas of CDM because of the generalized Liouville’s 
theorem that accounts for the time evolution of the metric. 
Under the influence of gravity, a particular phase space vol¬ 
ume of the hypersurface is stretched and folded with each or¬ 
bit of the CDM creating layers of fine-grained dark matter 
streams, each with a vanishingly small velocity dispersion. 
These stretches and folds also produced caustics (Tremaine 
1999), regions with very high CDM densities that are in¬ 
versely proportional to the square root of the velocity disper¬ 
sion, potentially providing large boosts to CDM annihilation 
into photons. Identifying such boosts is important because 
it would help explain the positron excess seen by PAMELA 
for energies above 10 GeV (Adriani et al. 2013), if such ex¬ 
cess stems from WIMP annihilation. Indeed, supersymmet¬ 
ric radiative corrections to the neutralino annihilation positron 
branching ratio can explain the rise in positron excess for ris¬ 
ing energies provided features (like caustics) in the halo struc¬ 
ture boost the dark matter density (Bergstrom et al. 2008). 

In order to understand how caustics could provide such 
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a boost, Vogelsberger and White (Vogelsberger and White 
2011) (henceforth V&W) integrated the geodesic deviation 
equation (Vogelsberger et al. 2008) in tandem with the equa¬ 
tions of motion of N simulation particles to analyze the 
fine-grained stream structure of the halos starting from gen¬ 
eral ACDM initial conditions. They were able to count all 
the caustics encountered by a simulation particle along its 
geodesic and calculate the CDM annihilation rate enhance¬ 
ments; they found that caustics provide at best a boost of 
0.1%. The V&W simulations also produced a probability 
distribution that a fine-grained stream containing a particu¬ 
lar fraction of the average local density would pass through a 
detector on Earth. For example, they found that there was a 
20% chance of a local CDM detection where a single stream 
accounted for 1% of the total signal. These results led V&W 
to conclude that the CDM velocity distribution was locally 
smooth and that the fine-grained halo structure would be too 
difficult to detect on Earth unless the CDM are axions. 

On the other hand, space-based detectors could discern the 
local fine-grained structure because the Earth acts as a grav¬ 
itational lens that separates dark matter streams 3 with differ¬ 
ent, sharply-peaked, velocities. In this paper, a Schwarzschild 
metric is used to show that weakly interacting particles 
streaming at 220 km/s (the approximate orbital velocity of 
the solar system around the galactic center (Karachentsev and 
Makarov 1996)) through a compact body with a small impact 
parameter (i.e., near the core) will experience huge CDM den¬ 
sity enhancements at the root (the nearest focal point com¬ 
posed of particles streaming through the core of the body) 
of the order of 10 9 for the Earth and 10 11 for Jupiter. For 
the Sun, the entire hair from root to tip (the focal point of 
particles streaming at the surface of a body ) 4 * * * is inside the 
solar radius for median CDM velocities of 220 km/s so that 
any direct detection near the Sun would be suppressed by the 
low probability of finding a high velocity dark matter stream, 
not to mention the challenges of close proximity to the Sun. 
Dense dark matter hairs provide an entirely unique opportu¬ 
nity to study dark matter interactions and local stream proper¬ 
ties. These CDM density enhancements will boost their direct 
detection rate for a space-bound experiment without invoking 
any physics beyond the standard model of particle physics, 
ACDM and general relativity and may be the only way to ob¬ 
tain empirical data on local dark matter streams which would 
further our understanding of halo structure and formation. 

This paper is organized as follows: In section 2, known re¬ 
sults for exterior geodesic solutions near Earth of both mass¬ 
less and massive particles are briefly reviewed for complete¬ 
ness and a description of the flux magnification for external 
solutions is given. Section 3 derives the analytic, constant 
density interior solutions of the geodesic equation (GE) and 
the density enhancement for massive particles is shown to be 
huge; the massive particle solution is expanded in powers of 
the impact parameter b , and the focal point is found to be 
independent of b to leading order with corrections of order 
(Z?/Re ) 2 where Re is the radius of the Earth (but could be 
the radius of any compact body). It is also shown that the 
non-zero but tiny velocity dispersion of the CDM has a neg¬ 

3 From hereon, dark matter stream refers exclusively to the primordial dark 
matter streams with tiny velocity dispersions, and not to tidal dark matter 
streams unless explicitly noted. 

4 It will be seen below that the flux of massive particles is also magnified 

for particles external to a body suggesting that the hair length may be infinite. 

However, the focal point distance increases too fast for external solutions 

undermining their experimental usefulness. 


ligible impact on the density enhancements of the dark mat¬ 
ter hairs for standard mass assignments of the CDM. In sec¬ 
tion 4, the Fast Accurate Integrand Renormalization (FAIR) 
algorithm for solving the GE for realistic radial densities is 
described and compared with the analytic constant density re¬ 
sult. Because FAIR is an algorithm that modifies the short 
scale radial dependence of the density profile, it can be used 
to demonstrate the robustness of the density enhancement re¬ 
sults. In section 5, the numerical solutions of the GE for a 
model Earth with a radial density given by the Preliminary 
Reference Earth Model (PREM) (Dziewonski and Anderson 
1981) are calculated using FAIR, and the hair focal points and 
density enhancements are plotted as a function of the impact 
parameter; the same plots are generated for Jupiter with a ra¬ 
dial density given by the model J1 l-4a from Ref. Nettelmann 
et al. (2012). In the last section, the implications of our re¬ 
sults are discussed, including the benefits of exploring a hair 
(which includes a powerful new tool to explore the interiors 
of compact bodies) and the conditions constraining a space 
mission dedicated to finding such a hair; an equation relating 
the different parameters of such a mission is also derived. 

2 . EXTERIOR SOLUTIONS AND DENSITY 
ENHANCEMENT FOR EARTH 

As a weakly interacting particle streams above the surface 
of the Earth, its trajectory curves towards a point on the axis 
going through the center of mass of the Earth and parallel to 
the incident velocity at t —> -oo. The location of this point 
and the corresponding magnification of flux intensity can be 
extracted from the known geodesic solutions for massless and 
massive particles in a Schwarzschild metric. The geodesic 
equation for a massive particle is given by 



while the equation for massless particles is 



where r s = 2GMe/c 2 is the Schwarzschild radius, Me is the 
Earth’s mass and G the gravitational constant. The r deriva¬ 
tive of 6 is zero thanks to angular momentum conservation 
with sin 6 = 1 . 

The exterior solution for a massless particle with impact 
parameter b > Re can be obtained by integrating 

^ =-b(l-b 2 u 2 + r s b 2 u 3 Y k (3) 

d u v ’ 

giving the usual solution 


0 = 0 ' + acos (bu) + 7 ^ , (4) 

b 



where u = 1 /r; this corresponds to a deflection angle 

** = 2 t 


( 6 ) 
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Fig. 1.— Simplified drawing (not to scale) of DM particles, belonging to a 
single stream, passing through a body. The hair is depicted in red, from the 
root (where the particles streaming near the core focus) to the tip (where the 
particles grazing the surface of a body of radius R focus). Figure 2 shows an 
inset of the square above, where dark matter particles with impact parameter 
frmin > b > b max focus near F(b). All the particles belonging to the same 
stream, and passing through the ring with perimeters located at b mm and b max , 
will focus on a circular detector with radius R&. 

Since r s = 9 x 10 -3 m for the Earth and b > Re = 6.371 x 10 6 
for geodesics external to the Earth, the deflection of mass¬ 
less non-interacting particles is effectively negligible with fo¬ 
cal points 


u 2 

F{b) = — > 2 x 10 15 m (7) 

2r s 

which are completely unreachable. The exterior solutions for 
massless particles with small impact parameters (b < Re) ex¬ 
terior to the Earth will be needed below for the initial angle 
value of the interior solutions of massless particles. They can 
be obtained by expanding the denominator in Eq. (3) and in¬ 
tegrating yielding 

<f> = n - ^ (l - Vl - b 2 u 2 ^ - arcsin(Z?w). (8) 

For massive particles neglecting the insignificant terms pro¬ 
portional to r s /r the exterior solution is the usual orbit equa¬ 
tion 
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To calculate the density enhancements by a planetary body 
of weakly interacting particles belonging to a dark matter 
stream 5 , a circular-disk detector with radius Re is taken to be 
located on the axis going through the center of the planet and 
parallel to v (see Figs. 1 and 2). In the Schwarzschild met¬ 
ric, every particle belonging to that dark matter stream will 
eventually cross that axis. The density enhancement M at a 
detector situated at point F(b) along the axis is given by the 
ratio of the annulus surface to the surface of the detector 

M[F(b)] = bmax ~ bmm = ^ . (14) 

5 All dark matter streams in this paper are defined by their velocity, v, 
and density ratio, p s /(p), where p s is the stream density and (p) is the local 
average density. 



F(b) 


Fig. 2.— All particles from a dark matter stream with impact parameter 
between b max and b m - m will cross the axis between F(b m \ n ) and F(b max ) and 
stream through (or hit) the circular detector of radius Ri ). 


Integrating over all possible contributions to a detector located 
at a point b yields 
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where the step function 6{x) = 0,1 for x < 0 and x > 0 respec¬ 
tively. Particle moving with speed v with b > Re experience a 
flux enhancement: 
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Although this is a large enhancement for b > Re, the focal 
point distance increases quadratically with the impact param¬ 
eter, making it very difficult to benefit from the linearly in¬ 
creasingly magnified flux of streaming massive particles. In 
addition, the interior solutions of the geodesic equation of 
massive particles streaming through the core will be found to 
be orders of magnitude larger, with focal points located within 
a million kilometers of the Earth and less for Jupiter. 


3. INTERIOR SOLUTIONS: CONSTANT DENSITY 
EARTH 


The interior geodesics of massless particles in a general 
Schwarzschild metric are solutions of 
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while the massive particles are solutions of 

i,2 \ „ i,2 u2 
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where ® is the interior potential of a body. For compact bod¬ 
ies with constant density, e~® ,c has an analytical form (Wald 
1984) 
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with p the average density of the body. 


Massless particle: Constant Density Interior Solution 



MASSLESS PARTICLES. In this analytical case, the 
differential equation for massless particles is 


d u 
d(p 
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with a solution that can be readily written down 
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Note that as b —> 7? E , 0' —> 0', and that as b —> 0, 0 { tt/2, 

as should be the case. Although the flux of massless 
particles streaming through a constant density Earth is greatly 
magnified as seen in Fig. 3, their focal points are located 
far outside the solar system (Fig. 4) making explicit the fact 
that planetary lensing and density enhancement of massless 
weakly interacting particles cannot be exploited locally 
for detection purposes 6 . In Fig. 5, the log-log plot of the 
enhancement as a function of the focal point distance is also 
given where the enhancement peak for massless particles 
streaming near the core is now clearly visible (it is pressed 
against the ordinate in Fig. 3). In addition, the approximate 
straight line in the log-log plot between the enhancement 
and the focal distance shows that any exponential density 
enhancement is accompanied by exponentially increasing 
distances. 


6 For massless or high energy neutrinos streaming through the Sun, the 
focal points are near the orbit of Uranus (Gerver 1988; Escribano et al. 2001) 


Fig. 3. — Density enhancement of massless particles streaming through a 
constant density Earth. Note that there is a barely visible peak for tiny impact 
parameter that is more salient in the log-log plot of Fig. 5. 


Massless particle: Constant Density Interior Solution 



Fig. 4.— Focal points as a function of impact parameter of massless particles 
streaming through a constant density Earth. 


MASSIVE PARTICLE. Using Eq. (20), Eq. (22) and 
Eq. (23), the radial equation of motion for a massive particle 
is given by 
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This is almost identical to Eq. (24) except for the c 2 /v 2 which 
now dominates and enhances the r s terms. The massive con- 
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Massless particle: Constant Density Interior Solution 



Fig. 5. — Density enhancement of massless particles streaming through a 
constant density Earth as a function of focal point distance. 

stant density interior solution is therefore 
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The corresponding density enhancement and focal point are 
plotted against impact parameter in Figs. 6 and 7 respectively. 
The massive particles streaming near the center of the planet 
clearly experience the largest magnification. To understand 
why, consider the expansion of (f)' M in powers of b/R^ 


, ^7T 3 b 3 b 3 

(Pm ~2 2 SR e R e + 8 ^ 4 e ' 

(34) 

c 3 b f, lb 2 ) 
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The critical thing to note is that the deflection angle 6(f) is lin¬ 
ear in b with corrections suppressed by factors of ( b/R e) 2 . As 
such, near the Earth’s center, 


which is consistent with Fig. 6 and independent of the DM ve¬ 
locity. Velocity independence is important because in a search 
for hairs, the density enhancement will only depend on the 
density ratio of the dark matter stream, p s /(p). In addition, 
the location of the peak density enhancement (which is about 
twice the value in Eq. (38)) is simply 

Vak = (4R d ReP = 55 km for /? D =1 m. (39) 

Note also that the singularity at Rd = 0 does not exist in real¬ 
ity because the non-zero velocity dispersion of the dark matter 
stream is a hard limit on the focus intensity. For a WIMP dis¬ 
persion of 3 x 10 -4 m/s assumed to be perpendicular to the 
direction of the WIMP, you get a total perpendicular distance 
travelled 

SRj) ~ d ± = —Fib) = 6 m at b ~ 10 km (40) 
v 

for a constant density Earth in the worst case scenario of an 
entirely perpendicular dispersion; for a 10 peV axion, d ± is a 
fraction of a millimeter. In other words, reducing the detector 
radius beyond SRv cannot result in a larger density enhance¬ 
ment and Eq. (38) can be rewritten 


M( 0) = 


mpRjp 

Rd + SRj) 


(41) 


The velocity dispersion at present age for WIMPs and axions 
is of order (Sikivie 1999) 

Sv a ~ 10“ 8 ( 10 5eV ) m/s (42) 

\ m a ) 

n I GeV 

£v w ~ 3 x 1(T 3 J -m/s (43) 

V m w 

It is therefore seen that for m a > 10 -9 eV and mw > 10 2 GeV, 
the density enhancement observed by a detector with a radius 
of a few meters will be little affected by the velocity disper¬ 
sion of the CDM. 

The fact that the deflection angle is linear in b/R e with cor¬ 
rections of order (b/R e) 3 is not unique to constant density pro¬ 
files as will be verified below for the Earth and Jupiter using 
realistic density profiles from the literature. We now turn to 
an algorithm that will permit us to calculate density enhance¬ 
ments for arbitrary density profiles. 


£R 2 ( u2 \ 

F(b) = 1 + —7 =l-6x 10 9 m, (36) 

3 1 4 Rl) 

which is consistent with the exact solution in Fig. 7. In other 
words, the focal point is effectively independent of the impact 
parameter b for particles streaming through the core within a 
radius of dozens of km’s resulting in a huge density enhance¬ 
ment. This fact is visually clear in Fig. 8 where the enhance¬ 
ment is plotted as a function of the focal point where the root 
of the hair would appear point-like to a space probe. To cal¬ 
culate the density enhancement near the core, it is only neces¬ 
sary to solve for b 

6cf>[Fib) - Fib)] = R d . (37) 

For example, taking b - 0 yields 


M(0) = 


(4 «)3 


Rd 


= 3 x 10 9 (R d = 1 m), (38) 


4. FAST ACCURATE INTEGRAND 
RENORMALIZATION (FAIR) ALGORITHM 

In order to calculate the geodesics for realistic planetary 
radial density models, it is necessary to quickly and accurately 
perform integrals of the form 


X 0f rn 

d + = 

i D/?body 


vb/(cr 2 )dr 
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for a general potential ® and where 0i is the initial angle de¬ 
termined by the exterior solution of the geodesic equation and 
where rf is the radius of closest approach at 0f. Since the 
integrand is singular at r - rf, 0f is difficult to evaluate nu¬ 
merically and the deflection angle may be hard to pin down to 
sufficient precision. 

Early in this work, Runge-Kutta methods were used to solve 
the geodesic equations, but they were found to converge too 
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Massive particle: Const. Dens. Interior Solution 



Fig. 6. — Density enhancements of massive particles streaming through a 
constant density Earth as a function of impact parameter b, assuming R d = 
1. Note the large spike at small impact parameters where the focal points 
are effectively independent of the impact parameter with corrections of order 
b 2 /R The density enhancement peak for a constant density planet can be 
located: for this Earth model, it occurs at impact parameter b - 55 km. 


Massive particle: Const. Dens. Interior Solution (V=220km/s) 



Fig. 7. — Focal points of massive particles streaming through a constant 
density Earth as a function of impact parameter b. The root of a constant 
density hair stemming from particles moving at 220 km/s is at a distance 
similar to L2 (with a different direction) and is relatively accessible. 

slowly. Other numerical methods found in the literature (Press 
2007) were tried with unsatisfactory results because 0f could 
not be calculated accurately enough near the caustic of the fo¬ 
cal point. Indeed, the density enhancement depends on the 
derivative of the impact parameter with respect to the focal 
point (see Eq. (15)) which is huge for small impact param¬ 
eters. To get around this problem, a method was developed 
to perform the integral in a quasi-analytic fashion that relied 
on the physical intuition that the precise details of the density 
should have little effect on the overall shape of the geodesics. 
In other words, the precise form of the density on a scale of 
a few meters should not significantly impact the interior solu¬ 
tion of a large body with a radius of thousands of kilometers. 
Additionally, none of the planetary models claim to be ac¬ 


Massive particle: Const. Dens. Interior Solution (V=220km/s) 



Fig. 8. — Density enhancements of massive particles streaming through a 
constant density Earth as a function of focal point distance. The enhancement 
curve is effectively vertical at the root of the hair suggesting a dense, point¬ 
like structure trailed by a sparser tail. Note that this ’’sparse tail” is still highly 
dense as compared to the ambient density of CDM. 


curate at small scales, which leaves us with the flexibility to 
choose a convenient density function that can be solved ana¬ 
lytically on the small radial scale 7 . This algorithm also allows 
a test of the robustness of the qualitative features of the density 
enhancement, like the sharp peak at low impact parameters. 

Assuming the prior existence of a radial density model for 
a particular body, the basic FAIR algorithm works as follows: 

• Split the integrand into N segments of length much 
smaller than the radius of the body. 

• Create a density profile within each segment that is an¬ 
alytically solvable. Equate the average density of the 
density profile in that segment to the known, average, 
model value within that segment. 

• Starting with the known fa+ \ = 0; from the exterior 
geodesics equation, use the analytic solution of the in¬ 
tegral inside that segment to calculate the next value of 
the angle 0*. 

• Do this until you reach the segment that contains the 
value of r = rf at some layer j. At that point, calculate 
fa with 

fa = (pj+i ~ (analytic solution evaluated at r - ry+i) (45) 

The FAIR algorithm was verified numerically by compar¬ 
ing the analytic solution of the constant density Earth with 
a highly unphysical case where each layer had a short scale 
radial dependence of 1/r 4 

■^layers - 1 

P= Z 7 (46) 

1=0 

7 The word renormalization refers to this flexibility since renormalization 
in quantum field theory is effectively a modification of short scale physics. 
For example, propagator functions can be modified by multiplicative func¬ 
tions that make integrands analytic at their short scale singularity. 
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Massive particle: Renormalization 1/r 4 Massive particle: Renormalization 1/r 4 residual 




Fig. 9.— Density enhancement results from using a r 4 radial dependence 
within the density segments. The black analytic solution is covered by the 
blue line of the FAIR algorithm with 10 7 layers. 


subject to the average-density-equality condition 
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where p/ is the density value extracted from a standard model 
at the / th layer. In the constant density case, pi = p for all 
V s. After a bit of algebra, the solution for each layer with the 
density profile of Eq. (46) is found to be 
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where w* = 1/r/ and 0/+i was calculated in the previous step. 
Comparing the results for the density enhancement of the an¬ 
alytic solution to the unphysical r -4 in Fig. 9, the agreement 
improved with the number of layers and was quite good for 
10 7 layers demonstrating a resolving robustness of the quali¬ 
tative rapid magnification rise for decreasing impact parame¬ 
ter. The residuals for both the density enhancement and focal 
point locations are provided in Fig. 10 further showing how 
the FAIR algorithm tends to the right answer with increasing 
number of layers even for a highly unphysical case. The 10 7 
layer case took about 10 minutes to run and demonstrates that 
even if your initial guess for the intra-layer density radial de¬ 
pendence is very off, you can still obtain a very good result 
compared to the true solution rather quickly. 

The simulated results appearing below for PREM and Jll- 
4a were obtained with the geodesolver code based on the 
above algorithm but with constant intra-density layers more 
in line with physical expectations. The iteration equations of 


Fig. 10.— Density enhancement (solid) and focal point (dash) residuals be¬ 
tween the constant density analytic and interior solutions solved using a r~ 4 
density renormalization of a constant density Earth for increasing number 
of layers. Using 10 7 steps, the r~ 4 solution is seen to be accurate to about 
b = 500 km which is orders of magnitudes better than the RK solution (not 
shown) for the same number of steps. 


geodesolver are 

dw\ 2 1 r 2 3 4l 

— = — \ao + aiu + a^u - u 

d (p) u z 1 J 

47tG 

a ° = ~3vW Pi 

+ ^2(l y '>ll + 47rG A' , il) 

2G (M^ 1 An _ 9 

a2 = w[—-Y piri 
Mi-1 = Mi - ~^Pi ( r , 3 - ff-i) 

Vi +1 = V i+ 2 + 4nGp i+ \ (rf +2 - rf +l ) 


(52) 

(53) 

(54) 

(55) 

(56) 

(57) 


where the starting values of r, M, and V are the radius of the 
planet, the mass of the planet, and zero, respectively; the p’s 
are given by the physical planetary model. To solve Eq. (52), 
it is best to remove the r 3 term with the following trick: 

u 3 = ul + Su 3 (58) 

= ul + 3 u 2 (u - uo) (59) 

and solve for u 3 . Eq. (52) now has the same form as Eq. (29) 
and can be solved in the same manner. 

The agreement between geodesolver and the analytic so¬ 
lution for the constant density Earth is excellent: using only 
3 x 10 5 layers, the density enhancements agree with the con¬ 
stant density analytical case at the 0.1% (Fig. 11) level while 
the focal points match to 6 significant digits (Fig. 12), includ¬ 
ing near the critical sharp density enhancement peak at small 
impact parameter. At 10 5 layers, running geodesolver takes 
only a few seconds to calculate each point. 


5. INTERIOR SOLUTIONS: PRELIMINARY 
REFERENCE EARTH MODEL AND J11-4A JUPITER 
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Massive particle: Const. Dens. Magnification Residual PREM Density 




Fig. 11.— Density enhancement residual between the analytic and 
geodes lover interior solutions of a constant density Earth for 3 x 10 5 layers. 


Fig. 13.— Preliminary Reference Earth Model radial density. 


Massive particle: Const. Dens. Focal Point Residual 



Massive particle: PREM Interior Solution 



Fig. 12. — Focal point residual between the analytic and geodeslover interior 
solutions of a constant density Earth. 


Fig. 14.— Density enhancement of massive particles streaming through a 
PREM model Earth as a function of their impact parameter. Earth’s layers 
are clearly visible showing that a hair could be a highly accurate probe of 
Earth’s interior. 


geodesolver was used on realistic planetary models for 
Earth (PREM (Dziewonski and Anderson 1981)) and Jupiter 
(Jll-4a (Nettelmann et al. 2012)); these density profiles are 
plotted in Figs. 13 and 18 respectively. The density enhance¬ 
ments and focal points for PREM are plotted in Figs 14, 15, 
and 16 while Figs 19, 20, and 21 are the corresponding plots 
for Jll-4a. The density enhancement plots are particularly 
striking as they clearly resolve the interior structure of the 
planets. Additionally, they also exhibit huge density enhance¬ 
ments for small Z?’s with point-like root structures as in the 
constant density Earth model. Other characteristics of the 
plots are tabulated in Table 1 . The second column of Table 1 
refers to the cutoff velocity of the hair, namely the velocity at 
which the root of the hair is below the surface of the planet. 
This velocity constrains searches for hairs surrounding mas¬ 
sive bodies as entire hairs from root to tip may find them¬ 
selves inside a body for likely values of the velocity (e.g., the 
Sun). That cutoff velocity can numerically be calculated us¬ 


ing geodesolver but it is also easily obtained from the roots 
already calculated at v = 220 km/s using the constant den¬ 
sity result that the root distance grows quadratically with the 
velocity (Eq. (36)): 



The fact that relations derived for a constant density Earth are 
applicable to realistic planetary models is explained below. 

To better understand the behavior of the focal 
points as b —> 0, let’s blow up the small impact parameter 
segment of the focal point plot in Fig. 15 and plot the deriva¬ 
tive [F(b i+ 1 ) - F(bi)]/(b i+ 1 - bi) as a function of b in Fig. 17. 
Thus, even for a realistic model such as the PREM, the deriva¬ 
tive of the focal point is linear with the impact parameter as 
b —> 0. This implies that the focal points are constants save 
for tiny corrections of the order Z? 2 /7^ ody . The fact that this 
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Massive particle: PREM Interior Solution (V=220km/s) 



Fig. 15.— Focal points of massive particles streaming through a PREM 
model Earth as a function of impact parameter b. The root of the hair is 
located at a distance of 165 Re- 


Derivative Of Focal Point Near Core Center 



Fig. 17.— Derivative of the focal point with respect to the impact parameter 
b for constant density (black curve) and PREM Earth (in red). Near b = 0, 
the derivative is dominated by numerical noise. 


Massive particle: PREM Interior Solution (V=220km/s) 



Fig. 16.— Density enhancements of massive particles streaming through a 
PREM model Earth as a function of focal point distance. 


J11-4a Jupiter Density 



Fig. 18.— J1 l-4a Jupiter radial density. 


result generally holds for arbitrary radial density profiles can 
be understood by rewriting Eq. (20) 


d 0 = 



vbdr 


rs M(r)]°- 5 U I O 

r M E J L ie 2 




( 61 ) 


and expanding in powers of b. Note that the singularity of the 
integrand is not within the integration interval, but at the lower 
limit where the integrated result is zero since this is the turning 
point where the radial distance is minimized, and about which 
the particle trajectory is symmetric. Thus, similar huge den¬ 
sity enhancements of coherent particles exiting the core can 
be expected for any compact body in the Schwarzschild met¬ 
ric, as was verified analytically for the constant density Earth 
and numerically for the PREM and Jll-4a. Quite generally, 
the focal points corresponding to impact parameters near the 
core of a Schwarzschild compact body can be parameterized 


by 


Fib) = F r + ^-Ir (62) 

In Table 1, these parameters are given in columns 3 and 4 for 
PREM and Jll-4a. 


6. DISCUSSION 

Key aspects of direct dark matter detection experiments that 
are not empirically constrained is the incident dark matter flux 
on Earth and the resulting lack of knowledge of the veloc¬ 
ity distribution of the CDM. Traditionally, the gravitation¬ 
ally bound CDM has been modeled as an ideal gas with a 
Maxwellian velocity distribution 

j exp 


/cdm(v) - 


AN 


(— 

yfnv S0 1 \v s0 l 


V 

V sol 


V < v e; 


(63) 
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Massive Particle: J11-4a interior solution 



Fig. 19.— Density enhancement of massive particles as a function of their 
impact parameter through Jupiter. Once again, the interior structure is clearly 
visible in the hair density enhancement profile. 


Massive Particle: J11-4a interior solution (V=220km/s) 



Fig. 20.— Focal points of massive particles streaming at 220 km/s through 
a J1 l-4a Jupiter model as a function of impact parameter b. The root of the 
hair is located at 3 Jupiter radius and any search seeking to take advantage of 
the large density enhancements of dark matter fluxes near Jupiter would also 
have to deal with the intense fields surrounding the gas giant. 


Body 

v c (km/s) 

FA m) 

U( m 1 ) 

PM 

ECD 

14 

mil 3 

f/6 


(4*0*1) 3 

*D 

T~ 

PREM 

18 

10 9 

4 x 1(T 5 

3.1 x 10 9 

Jll-4a 

125 

2.15 x 10 8 

1.3 x 1CT 6 

6.6 x 10 11 


TABLE 1 

Values for different planetary models: ECD are the values for an Earth 

OF CONSTANT DENSITY, PREM IS THE PRELIMINARY REFERENCE EARTH MODEL AND 

J1 1-4a is a realistic model for Jupiter. v c is the cutoff velocity below 

WHICH THE ROOT OF THE HAIR IS BELOW THE SURFACE OF THE PLANET. F R IS THE 
FOCAL POINT LOCATION OF THE ROOT OF THE HAIR FOR V = 220 KM/s, WHILE 
^(M" 1 ) PARAMETRIZES ITS b 2 DEPENDENCE (SEE TEXT). PM IS THE PEAK DENSITY 
ENHANCEMENT AT THE ROOT FOR EACH CASE. FINALLY, NOTE THAT FOR PREM, THE 
GREATEST DENSITY ENHANCEMENT OCCURS FOR & PEAK =48 KM, SIMILAR TO THE 
RESULT ANALYTICALLY DERIVED FOR THE ECD MODEL (Eq. (39)) AND 
^peak- 1133 km for Jupiter. 


Massive Particle: J11-4a interior solution (V=220km/s) 



Fig. 21.— Density enhancements of massive particles streaming at 220 km/s 
through a J1 l-4a Jupiter model as a function of focal point distance. 

where /V is a normalization constant to ensure that the integral 
over the distribution results in the local CDM density, v so i is 
the Sun’s velocity with respect to a CDM halo taken to be sta¬ 
tionary, and v esc is the escape velocity of the CDM from the 
galaxy. This isotropic, analytically smooth velocity distribu¬ 
tion is unrealistic since the CDM had tiny primordial velocity 
dispersion at the last scattering surface. It only acquired an¬ 
gular momentum and velocity dispersion through their inter¬ 
actions with gravitational wells and tidal forces. The fraction 
of CDM that would have been thermalized and acquired an 
approximate Maxwellian velocity distributions are those that 
would have gone through many orbits about the Milky Way. 
There will also be a fraction of CDM particles that has un¬ 
dergone fewer orbits such that the velocity distributions has 
peaks (Sikivie et al 1995). 

A number of realistic theoretical and numerical calcula¬ 
tions have been performed to try and understand the veloc¬ 
ity distribution in the halo (Springel et al 2008; Hansen 
2009; Vogelsberger et al. 2009). In particular, comprehen¬ 
sive and realistic fine-grained simulations (Vogelsberger and 
White 2011; Springel et al 2008; Vogelsberger et al 2009) 
of the geodesics followed by the CDM produce a distribution 
of finely parsed streams with discrete velocities that can be 
treated effectively as smooth for experiments on Earth’s sur¬ 
face. There are millions of streams flowing through our solar 
system with p s > 10 _7 (p), including at least a few streams 
with at least 10 _3 (p) (Vogelsberger and White 2011). De¬ 
tecting hairs stemming from these streams would be a huge 
scientific windfall as they would offer us 

• a nearby, extraordinarily intense source of cold dark 
matter particles accessible to space detectors studying 
the nature and interactions of CDM; 

• our only way to directly explore the local fine-grained 
stream structure, a prediction of ACDM; 

• a direct and detailed look at the radial density distribu¬ 
tion of any planet or moon, since determining the loca¬ 
tion of a single hair near Earth or Jupiter would imme¬ 
diately provide the location of that same hair extending 
from any other body in the solar system. The densest, 
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fine-grained dark matter streams, that have undergone 
the fewest orbits, can be expected to be very smooth on 
length scales similar to the solar system. Indeed, the 
original, unclustered, dark matter flows that stretched 
and folded into fine-grained dark matter streams, would 
have been much larger than a solar system Natarajan 
and Sikivie (2008). The original dark matter stream 
would have had a CDM density near the cosmological 
average, but the subsequent stretching and folding ex¬ 
plains why the particle density in a fine-grained stream 
is typically small: the total number of collisionless dark 
matter particles remained constant in an expanding vol¬ 
ume. The densest dark matter streams should therefore 
have the same velocity and density within the solar sys¬ 
tem resulting in hairs with near identical orientations at 
different planets 8 . 

Locating these hairs is likely to be challenging but may be fa¬ 
cilitated by the fact that the density enhancement is indepen¬ 
dent of the stream velocity, with minimum density enhance¬ 
ments of order 10 7 for Earth and 10 8 for Jupiter. In addi¬ 
tion, the detailed spherically symmetric density structure of 
the planetary layers does not impact the qualitative features 
of a hair such as a sharp rise at low impact parameters. 

Although planets are not precisely described by a 
Schwarzschild metric, slight deformations should not negate 
the basic arguments presented in this paper. The reason is that 
planetary deformations due to rotational motion do not break 
the mirror symmetry of the body across its center of gravity. 
This mirror symmetry constrains the metric to depend on the 
square of the angular momentum, since a particle with an¬ 
gular momentum L and another with angular momentum -L 
will meet at the same focal point. Hence, although planets 
are not precisely spherically symmetric, the hairs will reflect 
those deformations in a proportional fashion such that a hair 
cross-section will simply represent a 2-dimensional image of 
the planet from that vantage point. The issue of planetary de¬ 
formation will be explored numerically in a future paper with 
a generalization of the FAIR algorithm. 

Turning now to the parameters necessary for an estimate of 
the time and resources required for finding a hair, we have: 

• p represents the percentage of the sphere required for 
the search area. This number needs to be calculated 
since the orientation of the hairs is biased by the galac¬ 
tic orbital motion of the solar system. In addition, hairs 
with relatively large densities (say p s /(p ) > 10 -4 ) have 
likely gone through fewer orbits about the galaxy, and 
may have a different orientation distribution from the 
low density streams. 

• ftprobe is the number of probes simultaneously conduct¬ 
ing a search. In estimating the time required to find a 
hair, ft pro be will be solved for as a indication of the fea¬ 
sibility of the search. 

• f t is the average transit time of a probe through a hair. 
This parameter is critical for the strength of the signal 
in the event the probe does fly through a hair. If the 
transit time is too short for the sensitivity of the probe, 
the hair will not be detected. 

• h w , the average hair width crossed by a transiting probe. 

8 Any observed differences in the hairs between different planets would 
stem from local corrections of the Schwarzschild metric due to nearby bodies. 


• dprobe the distance from the body where the probe is 
concentrating its search. In practice, one would expect 
dprobe - to afford the best odds. 

• N s the number of streams for a particular density ratio 
Ps/<P>- 

The parameters above can be used to estimate the search time. 
To solve for ft pro be> that search time must be limited in some 
fashion. If one requires the search time to be less than a few 
years, the following equation follows 

rji 

search time = —— - _^ r ° be < 10 8 s (64) 

^probe^s h w 

£ = - ^ _ < 1CT 11 for Earth or (65) 

^probe N s h w m 

< 1CT 9 -At for Jupiter. (66) 

m z 

The parameter £ depends on the details of the experiment, like 
the choice of target, the target size and the expected event rate. 
In addition, although the transit time is constrained by the de¬ 
tector sensitivity, in practice it will also depend on the details 
of the probe orbit around the planetary body and how much 
fuel would be available to maneuver it around a desirable area 
of the focal sphere. These experimental parameters will in 
turn determine the density enhancement required to detect a 
hair during an average transit time, T t . Taking as a require¬ 
ment that a probe crossing a hair of width h w see the same 
number of detection events as an Earth-bound mission would 
see in a year for an equivalent experimental setup, a constraint 
on T t can be imposed 

7 ^-Mft ~ 10 7 s with p-M > o (67) 

<P) (P > 

where M is the average dark matter flux density enhancement 
experienced by the probe of radius R D as it transits across the 
width h w of a hair and o parametrizes the transit time and is 
therefore related to the probe perpendicular velocity through 
the hair cross-section. Noting from Fig. 10 of Ref. Vogels- 
berger and White (2011) that (pJ(p))N s ~ F(> p s ), £ can be 
rewritten 

p 10 7 s 

f- — - (68) 

Vobe^w MF{> Ps) 

where F(> p s ) is the probability that a dark matter particle 
belongs to a stream with a density greater than p s . As de¬ 
scribed above, solving this equation involves a number of un¬ 
knowns that are beyond the scope of the current paper which 
is dedicated to describing the new ACDM prediction of con¬ 
centrated dark matter hairs projecting from compact bodies. 
A comprehensive solution for ft pro be would look at a realistic 
experimental setup and provide likelihood parameter figures 
drawn from a sample of simulated planetary hair realizations 
as well as incorporate orbital/fuel constraints. However, one 
can get a rough idea of the difficulty of finding a hair by ap¬ 
proximating M ~ M/h w , and taking M ~ 10 8 for Jupiter and 
M ~ 10 7 for Earth to obtain 

ftprobe > — ° p 7 for Jupiter (69) 

p s F(> p s ) 

n probe > 10 4 ~~~~ r for Earth (70) 

Ps F(>p s ) 
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The reason p s appears in Eqs. (69,70) is because the number of 
streams at the smallest densities dwarfs the number of streams 
with the highest densities. For example, the simulations of 
Ref. Vogelsberger and White (2011) see 10 14 streams near the 
Sun but only about 10 6 are massive. The way to interpret the 
above equations is as a lower limit on the number of probes 
required to see hairs generated from streams with density p s . 
With this understanding, it is seen from Fig 10 of Ref. Vogels¬ 
berger and White (2011) that the ratio (p)/[p s F(> p s )] = 10 5 
for the densest streams (p s /(p) > 10 -4 ) with the ratio increas¬ 
ing very fast with decreasing p s . In particular, the least dense 
streams are effectively invisible because their hair width h w 
is tiny as seen from Eq. (67). For the dense streams with fi¬ 
nite /z w , ftprobe seems prohibitively high especially since the 
parameter o can be large. 

This rough estimate based on the simulations of Ref. Vo¬ 
gelsberger and White (2011) and assuming Earth-based de¬ 
tector sensitivity therefore appears to put hair detection out 
of reach, unless 1) p can be constrained to be quite small, 
2) the uncertainty on the radial distance to the hair root can 
be narrowed enough that density enhancements of 10 8 - 10 9 
for Earth and 10 10 - 10 11 for Jupiter can be used above, sup¬ 
pressing ftprobe by a factor of 10 2 - 10 4 , 3) future theoreti¬ 
cal/numerical developments allow the high-p s slope of Fig 10 
from Ref. Vogelsberger and White (2011) to be steeper 9 and 
4) a way is found to exploit the characteristics of the space en¬ 
vironment (such as microgravity) to design novel direct dark 
matter detectors, for both WIMPs and axions, with better sen¬ 
sitivity. 

One highly developed technology for space missions is re¬ 
frigeration for experiments requiring very low temperatures 
typical of axion searches (Asztalos et al. 2010). The cosmic 
microwave background (CMB) Planck satellite mission relied 
on active cooling to lower the temperature of the High Fre¬ 
quency Instrument bolometer plate to 93 mk (Ade et al. 2011) 
which is similar to the temperature goal for the ADMX appa¬ 
ratus upgrade. Generally speaking, axion searches and CMB 
space telescopes broadly share a reliance on microwave sci¬ 
ence and engineering which would make it more likely that 
a space-based axion detection mission could recycle well- 
tested CMB technology. Fow temperature detectors could 
also be efficaciously used for a WIMP detection space mis¬ 
sion thanks to their excellent signal-to-noise ratio and their 
event-by-event background rejection capabilities which may 
reduce the cosmic rays shielding requirements (Mirabolfathi 
2013; Guo and McKinsey 2013). Since discovering a hair is 
a separate goal with different benefits from direct dark matter 
detection (which may happen on Earth first), it may be possi¬ 
ble to transfer improvements in detector sensitivity to a space 
mission that would greatly increase the odds of finding a hair. 

Turning next to whether the compact body hairs could con¬ 
tribute a measurable signal to indirect dark matter searches 
from dark matter annihilation, consider the usual ambient 
contribution estimated from the equation 


0 (£’ r , if/) 


(crv) 

87rm 2 


S(E - E y ) 



p 2 (l , i//)dl , 


(71) 


where a monochromatic y-ray decay was assumed, m x is the 
CDM mass and the astrophysical content of the flux <D is con¬ 
tained in the line-of-sight integral on the end. Including the 


9 Specifically, Monte Carlos simulations should be run using the Planck 
cosmological parameters (Ade et al 2015), different gravitational softening 
values and narrower radial shells to estimate error bars on the slope. 


planetary focussing, the astrophysical integral becomes 

f p 2 (L >p)dl ^ f p 2 (/,^x// + G L (0,-A) (72) 

Jlos Jlos 

Cl (0, 0) = f d/Qcoft 0? 0) 

Jlos 

x X f dtdA(t)pl an .(t, Lip) , (73) 

hair Jhair 

where Elcoih 0,00 is the density of compact objects in the 
direction (0, 00 that focus the CDM beyond their surface, 
Phair(£, U 0) is the CDM density at point t along the hair of the 
gas giant, a distance / away from us with an angle 0 from the 
galactic center and /\air is the probability of seeing a hair with 
density phair; the last integral is performed about the volume of 
the hair with dA the area of the cross section at t. It is imme¬ 
diately clear why that second contribution is negligible: Qco 
has units of kpc -2 and can be completely ignored in CDM 
annihilation searches. This conclusion also holds locally as 
the number of annihilation products produced by each unit 
volume within the solar system far exceeds the number pro¬ 
duced along the hairs of planetary bodies even after account¬ 
ing for the density enhancement. This result also negates the 
possibility of detecting dark matter annihilation in hairs sur¬ 
rounding a particular compact body (such as a white dwarf or 
a neutron star) with a telescope, as it would be swamped by 
the ambient dark matter annihilation not to mention the direc¬ 
tional suppression stemming from the tiny opening angle of a 
created photon propagating from the compact body to an ob¬ 
serving satellite orbiting the Earth. This is also true for solar 
system bodies like Jupiter and Earth, which cannot produce 
detectable amounts of annihilation products. 

Lastly consider the possibility of relic neutrino hairs. The 
relic neutrino energies follow a Fermi-Dirac distribution and 
cannot form the caustics necessary to produce hairs, since 
their velocity distribution is completely smooth with a large 
dispersion. 

SUMMARY. A new ACDM prediction of concentrated 
dark matter hairs extending from compact bodies was 
described. The hairs are caustics of dark matter streams 
with primordial velocity dispersion passing through compact 
bodies. These caustics appear for both massless and massive 
particles but only the massive particles have caustics located 
near the lensing planet. It was shown that the existence of 
hairs is robust for realistic radial density profiles. It was also 
shown that the density enhancements are neither affected by 
the detailed structure of the body’s radial density profile, nor 
the dark matter stream velocity. Hair density enhancements 
were calculated for 3 cases: the analytic constant density 
Earth, the PREM and the Jll-4a density model of Jupiter. In 
all cases, large density enhancements were shown to exist up 
to 10 9 for Earth and 10 11 for Jupiter. Furthermore, the density 
enhancements of the hairs clearly reflected the boundary 
layers of the density profiles of the planets. The discovery 
of a hair would therefore be a huge boon for both planetary 
science and cosmology: 

• A hair would provide a uniquely powerful laboratory to 
study dark matter interactions. 

• Hairs are unique windows into the fine structure of the 
local dark matter streams. 
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• A hair discovered near Earth (or Jupiter) would find 
counterparts in all solar system bodies providing a pow¬ 
erful universal tool to probe the interior of almost any 
planet or moon. 


A set of parameters determining the feasibility of finding 
a hair was listed and an approximate equation for the 
likelihood of discovering a hair was written down, given 
certain constraints and assumptions. Implications for indirect 
dark matter searches and relic neutrinos were shown to be 
inexistent. 
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